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SECTION  I 


INTRODUCTION 


The  human  body  can  always  be  damaged  by  sufficiently  high  acceleration  levels.  Even  the 
gravitational  acceleration  of  32.2  ft/sec  (ref.  1)  can  cause  injuries  (damage  to  the  circulatory 
system,  for  example)  in  some  individuals,  after  a  sufficient  period  of  time  has  elapsed. 

The  higher  acceleration  levels  associated  with  re-entry  of  a  space  vehicle,  which  are 
maintained  for  only  a  few  minutes,  give  rise  to  "hydraulic  problems,"  caused  by  the  dis¬ 
placement  of  blood  and  other  fluids  from  their  normal  location  in  the  bod} .  In  contrast, 

"brief  acceleration"  or  "short  period"  acceleration  which  endures  for  a  period  of  less  than 
a  second,  can  cause  structural  damage,  if  the  level  is  high  enough,  but  cannot  significantly 
distort  fluid  flow  patterns.  It  is  this  short  period  category  of  accelerations,  such  as  those 
that  occur  in  aircraft  escape  systems,  and  the  body's  response  to  it,  which  is  the  subject 
of  this  report. 

The  precise  nature  of  the  structural  Injury  caused  by  short  period  acceleration  depends  up¬ 
on  the  orientation  of  the  acceleration  vector  with  respect  to  the  human  body,  the  restraint 
of  the  human  body,  the  magnitude  of  the  acceleration  and  the  way  in  which  the  magnitude 
varies  with  time.  This  last  aspect  is  sometimes  referrred  to  as  the  "frequency  content" 
of  the  input.acceleration,  but  such  terminology  can  be  misleading.  A  reference  to  "frequency 
content"  implies  that  such  descriptors  as  the  power  spectral  density  of  an  input  acceleration 
pulse  can  be  correlated  with  injury  mode,  and  at  the  present  time  this  is  not  proven;  indeed, 
quite  cogent  arguments  can  be  made  for  the  contention  that  power  spectral  density  is  not 
correctable  with  injury,  at  least  in  a  simple  way.  . 

The  basic  problem,  to  be  studied  can  be  summarized  as  follows.  Given  an  acceleration 
time  history,  such  as  that  associated  with  the  initial  trajectory  of  an  escape  system,  will 
it  cause  injury  when  applied  to  a  human  subject,  and  what  will  be  the  nature  of  the  injury  ? 

Of  course,  there  is  a  simple  and  direct  way  of  determining  the  answer  to  this  question  by 
experiment,  but  we  are  prevented  from  using  it  because  of  obvious  moral,  legal,  humani¬ 
tarian  and  other  prohibitions.  Yet  we  cannot  avoid  answering  the  question  because  we 
could  then  be  condemning  some  of  our  fighting  men  to  death  or  injury  in  the  future,  either 
because  they  will  have  to  use  escape  systems  with  excessively  dangerous  acceleration 
characteristics,  or  because  their  escape  envelope  is  unnecessarily  restricted  in  relation 
to  the  total  envelope  of  their  vehicle.  And  in  a  broader  context,  we  must  also  find  an 
answer  to  the  question  because  of  the  thousands  of  people  who  unnecessarily  lose  their 
lives  in  automobile  accidents  and  other  crash  environments;  deaths  which  might  often  be 
avoided,  given  adequate  restraint  and  protection. 

It  is  generally  agreed  that  the  best  line  of  attack  is  to  develop  an  accurate  description  of 
the  human  body,  including  both  its  dynamic  characteristics  and  the  static  load  bearing 
capability  of  its  many  complex  components,  so  that  we  can  accurately  calculate  the  effect 
of  a  given  input  acceleration  time  history,  and  also  that  we  may  be  able  to  design  a  re¬ 
straint  or  support  system  which  optimumly  couples  it  to  the  driving  force  associated  with 
the  acceleration.  The  construction  of  such  a  "dynamic  model"  is  a  task  routinely  under¬ 
taken  in  the  aerospace  industry  (among  others)  for  structures  which  are  built  from  well 
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understood  materials.  A  simple  example  is  an  aircraft,  for  which  the  strength  of  each 
element  is  accurately  calculated,  in  order  to  ensure  that  it  can  perform  its  function  with¬ 
out  failure.  We  also  determine  the  deformation  of  the  aircraft's  structure,  and  construct 
dynamic  models  which  tell  us  how  it  will  behave  in  response  to  gusts,  how  aeroelastic 
deformation  influences  its  stability  and  handling  characteristics,  and  to  assure  that  flutter 
and  divergence  cannot  occur  with  the  flight  envelope.  But  such  a  task  is  very  much  harder 
in  the  case  of  the  human  body,  firstly  because  it  is  so  much  more  complex  than  an  aircraft, 
and  secondly  because  many  of  its  constituent  materials  are  not  "simple"  linear  materials 
obeying  Hooke's  law  and  having  small  damping.  Additionally,  the  dynamic  characteristics 
can  be  varied  by  tensing  or  relaxing  the  muscles  which  make  up  so  much  of  the  body. 

Other  variations  occur  because  the  various  body  materials  tend  to  change  their  charac¬ 
teristics  with  age,  environment,  and  the  degree  of  physical  fitness  of  the  subject,  and  be¬ 
cause  people  are  not  manufactured  to  the  same  close  tolerances  as  aircraft. 

At  first  sight  one  might  suppose  that  an  adequate  description  of  the  human  body  could  be 
built  up  by  measuring  the  various  quantities  of  interest  from  cadaver  material.  Informa¬ 
tion  of  this  type  is  certainly  helpful,  and  indeed  the  pioneer  of  biodynamics,  Sigmund  Ruff 
(ref.  2),  obtained  a  remarkably  accurate  estimate  of  the  ability  of  the  human  spine  to  with¬ 
stand  ejection  seat  acceleration  in  this  way.  It  does  have  severe  limitations,  however, 
because  cadaver  material  is  hard  to  obtain  for  the  age  groups  of  primary  interest,  and  its 
dynamic  and  even  static  properties  are  thought  to  be  substantially  different  from  the  in  vivo 
case.  Additionally,  the  important  contributions  of  the  muscle  structure  and  fluids  are  al¬ 
most  entirely  absent  in  the  cadaver. 

It  is  therefore  necessary  to  devise  experiments  with  live  human  subjects  from  which  mean¬ 
ingful  data  can  be  obtained.  Since  it  is  essential  to  avoid  injury  to  the  subjects,  the  accel¬ 
eration  levels  used  in  such  tests  must  be  substantially  lower  than  those  which  would  cause 
injury,  and  the  results  have  to  be  extrapolated  in  some  way  in  order  to  predict  the  level  at 
which  injury  would  occur.  It  is  also  necessary  to  devise  tests  which  will  enable  us  tp  de¬ 
termine  the  damping  in  the  various  elements  of  the  system,  and  the  primary  dynamic  char¬ 
acteristics,  such  as  resonant  frequencies.  Orie  gross  way  of  doing  this  is  to  excite  the 
subject,  either  by  a  steady  state  sinusoidal  vibration,  or  by  an  impact,  and  to  measure  the 
relative  motions  of  various  parts  of  his  anatomy.  This  is  a  valuable  technique  in  the  posi¬ 
tive  spinal  case  (for  example)  where  relatively  large  and  rigid  masses,  such  as  the  head 
and  upper  torso,  are  moving  on  the  neck  and  spine  respectively.  Such  measurements  have 
been  attempted  both  optically  and  by  means  of  accelerometers  mounted  on  the  subject's  ' 
body.  Such  "external"  observations  cannot  tell  us  much  about  the  movements  of  the  various 
organs  inside  the  body,  of  course.  Alsc,  while  it  is  quite  easy  to  attach  a  transducer  to 
the  hip  or  the  upper  shoulder,  it  is  more  difficult  to  ensure  that  this  instrument  accurately 
tracks  the  motion  of  the  component  to  which  it  is  attached.  * 

One  particularly  simple  method  of  studying  the  app'areut  dynamic  response  of  the  human 
body  is  to  measure  its  driving  point  impedance  under  steady  state  sinusoidal  excitation. 
Because  no  instrumentation  is  attached  to  the  body,  this  technique  was  at  one  time  re¬ 
garded  as  nearly  ideal.  In  fact,  from  a  purely  experimental  point  of  view  it  is  very  ele¬ 
gant,  and  were  it  not  for  the  fact  that  impedance  measurements  taken  in  this  way  vary  by 


*  See  Goldman  and  Von  Gierke  in  reference  1,  for  example. 
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as  much  as  100%  from  subject  to  subject,  and  even  vary  for  a  given  subject  in  successive 
tests,  the  experimenter  would  be  well  satisfied  with  it.  But  such  data  turns  out  to  have 
rather  limited  value  to  the  theoretical  blodynamicist  in  his  construction  of  a  dynamic  model. 
As  will  be  shown  in  this  report,  impedance  data  is  of  little  value  because  it  is  only  meaning¬ 
ful  when  the  system  is  composed  of  linear  springs  and  dampers.  Additionally,  impedance 
turns  out  to  have  very  "limited  visibility"  in  a  series  coupled  system  with  finite  damping 
elements.  Putting  this  another  way,  we  might  regard  an  impedance  measurement  as  a 
"ray  of  light"  which  permits  us  to  'look  into"  the  human  body.  But  we  then  find  that  we 
~re  looking  into  a  rather  opaque  medium.  The  ray  of  light  shone  up  from  the  buttocks  (as 
in  the  measurement  of  driving  point  impedance  for  a  seated  subject)  tells  us  a  considerable 
amount  about  the  first-degree-of-freedom,  namely  the  pelvic  mass  moving  on  the  spring  of 
the  buttocks.  However,  it  gives  us  very  little  Illumination  into  the  dynamic  system  imme¬ 
diately  above  the  pelvic  girdle,  and  for  practical  purposes,  none  for  any  of  the  systems 
above  that. 

One  object  of  this  report  is  to  examine  critically  the  impedance  techniques  as  a  means  of 
obtaining  dynamic  information  about  the  human  body  by  simulating  the  same  measurements 
on  a  mathematical  model,  the  coefficients  of  which  are  known.  We  then  put  ourselves  in 
the  position  of  an  experimenter  who  is  trying  *o  determine  the  nature  of  the  model  from  his 
external  measurements  alone,  on  order  to  see  whether  he  would  meet  with  any  success.  In 
general,  we  conclude  that  he  would  not,  unless  all  the  dynamic  elements  in  the  body  were 
linear,  and  the  damping  was  very  low. 

If  there  is  an  overall  conclusion  to  this  work,  it  is  this;  there  is  no  "right  way"  to  discoyer  - 
the  "right  model"  for  the  human  body.  We  must  laboriously  build  a  model  piece- by-piece, 
sometimes  relying  upon  intuition,  sometimes  on  tht  biologist,  and  sometimes  relying  bn 
the  results  of  a  fortuitous  experimental  observation,  but  constantly  comparing  the  apparent 
external  behavior  of  our  model  with  experimental  observations  such  as  those  obtained  on 
the  AMRL  drop  test  facility  in  order  to  make  sure  that  it  is  not  straying  from  the  reality 
of  physical  measurements.  It  will  be  a  long  time  before  we  can  feel  that  we  have  a  com¬ 
pletely  adequate  model,  and  many  hundreds  of  workers  will  then  have  contributed  to  its 
construction,  just  as  in  any  other  discipline.  In  the  present  report  we  have  developed  a 
four-degree-of-freedom  model  of  a  seated  man  subjected  to  vertical  accelerations;  The 
characteristic  parameters  of  the  model  are  selected  by  making  use  of  relevant  available 
data  and  by  matching  calculated  driving  point  impedance  characteristics  with  measured 
values.  It  is  realized  that  a  number  of  multidegree-of-freedom  models  have  been  developed 
previously,  but  it  is  believed  that  development  of  the  model  in  this  way  allows  some  impor¬ 
tant  conclusions  to  be  reached  concerning  the  dynamics  of  the  body. 
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SECTION  II 


THE  CONCEPT  OF  MECHANICAL  IMPEDANCE 


In  structural  engineering  the  concept  of  the  linear  dynamic  system  is  widely  used.  Since 
many  structures  are  made  up  of  metal  components,  each  of  which  obeys  Hooke’s  law  (the 
linear  relationship  between  stress  and  strain)  the  dynamic  response  of  the  whole  is  quite 
accurately  described  by  a  matrix  of  second  order,  linear  differential  equations,  of  the 
type  used  to  describe  linear  electrical  circuits.  In  the  most  simple  example,  the  motion 
of  a  single-degree-of-freedom  mass/spring/dash  pot  mechanical  system  driven  by  a 
forcing  function  applied  to  the  mass  (see  figure  la)  can  be  described  by  the  equation 

my  +  2  Ky  +  ky  =  F(t) 

or  mv  +  2  Kv  +  kjfvdt  =  F(t) 


This  is  analogous  to  a  simple  electric  circuit  in  which  a  resistance,  a  condenser,  and  an 
inductance  are  connected  in  scries,  in  which  the  relationship  between  current  and  input, 
voltage  is  given  by: 

t  ‘  -  '' 

LI  +  HI  +  -=-/ldt  =  E(t)  ft 


Thus  for  any  consistent  sets  of  units  there  is  a  direct  correspondence  between  the  follow¬ 
ing  quantities:  '  * 


mass  m 

spring  constant  k 
damping  constant  ZK 
velocity  v=y 
driving  force  F(t) 
time  t 


Inductance  L  j 

reciprocal  of  capacitance 

resistance  R 

current  I 

input  voltage  E(t) 

time  t 


It  is  therefore  natural  to  take  advantage  of  the  powerful  analytical  techniques  developed 
by  electrical  engineers  over  the  last  eighty  years.  The  concept  of  impedance  is. one  of 
these ‘techniques,  which  for  the  simple  case  shown  in  figure  la  translates  into  mechan¬ 
ical  system  technology  as 

Electrical  Impedance  =  Z  =  R  +  i(£?L  + 

F 

and  Mechanical  Impedance  =  Z  =  =  2K  +  i(£m  -  k'fl) 

c 


so  that,  in  this  case,  the  real  and  complex  parts  of  Z  are  given  by: 

R(Z)  =  2K  and  I(Z)  =  Qm  -  k/£ 
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The  complex  notation  is  the  most  convenient  way  of  expressing  the  impedance  in  terms  of 
both  its  magnitude  and  phase  relationship.  Figure  1  shows  diagramatically  how  the  imped¬ 
ance  characteristics  of  a  number  of  mechanical  systems  will  depend  on  the  characteristics 
of  the  system.  It  should  be  observed  that  the  characteristics  of  those  parts  of  the  system 
closest  to  the  driving  point  are  fairly  readily  discernible  from  the  impedance  characteris¬ 
tic,  but  that  the  fuvther  from  the  driving  point  they  are  the  less  influence  they  have  on  the 
impedance  characteristic.  Mechanical  impedance  of  a  system  may  be  calculated  from 
direct  measurements  of  the  driving  force  magnitude  and  velocity.  Dynamically  similar 
systems  have  equal  impedance/frequency  characteristics  so  that,  inversely,  if  the  im¬ 
pedance  of  two  systems  are  equal  then  there  is  good  reason  to  suppose  that  they  are 
dynamically  similar. 

At  first  sight,  this  technique  seems  attractive  as  a  means  of  studying  the  dynamic  response 
of  the  human  body.  Alternative  methods,  such  as  attaching  accelerometers  to  various  parts 
of  the  test  subject's  anatomy  have  well  known  shortcomings;  in  particular,  it  is  difficult  to 
mount  them  "rigidly, "  so  that  they  move  with  the  skeletal  structure,  because  they  are 
necessarily  attached  to  the  relatively  soft  flesh  covering  the  structure.  It  is,  therefore, 
easy  to  understand  the  attraction  of  a  test  procedure  which  requires  no  instrumentation 
of  the  subject,  but  merely  an  unambiguous  measurement  of  input  force  and  velocity. 

Unfortunately,  as  indicated  in  figure  1  and  as  will  be  shown  later  in  this  report,  the  driving 
point  impedance  of  a  series  connected  lumped  parameter  model  only  gives  useful  informa¬ 
tion  for  those  few  systems  which  are  closest  to  the  driving  point. 

For  the  series  system  illustrated  in  figure  Id  for  example,  measurement  of  the  driving 
point  impedance  will  show  a  large  peak  corresponding  to  the  resonance  of  spring  and  damp¬ 
ing  k,K,  in  which  all  the  other  masses  tend  to  move  with  mass  m.  The  larger  the  damping 
ratios,  the  less  effect  systems  further  from  the  driving  point  have  on  the  impedance  charac¬ 
teristic  and  the  more  the  features  of  the  characteristics  will  be  obscured  (see  fig.  16). 

An  apparent  frequency  and  damping  can  be  deduced  for  this  mode,  but  it  bears  little  resem¬ 
blance  to  the  actual  values  of  system  (1),  or  Indeed  tc  any  simply  identifiable  characteristic 
of  the  model.  Thus,  the  procedure  sometimes  used  of  comparing  experimental  impedance 
measurements  with  a  s ingle- degree-of- freedom,  lumped  parameter  mode,  may  be  very 
misleading.  Much  more  sophisticated  methods  of  data  analysis  are  required  to  recover 
any  useful  information  from  impedance  data,  starting  with  computer  techniques  of  smooth¬ 
ing  the  experimental  data  in  a  mathematically  meaningful  way.  Even  when  this  is  done, 
however,  only  the  first  system  in  the  series  is  clearly  visible,  when  the  damping  is  as  high  ' 
as  that  of  the  human  body.  * 


A 


B 


C 


Mechanical  Equivalent  Mechanical  Magnitude  of 

System  Electrical  Impedance  Impedance 

System  Z  e  Force/Velocity 
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A  TWO-DEGREE-OF- FREEDOM  SYSTEM 

Figure  2  shows  the  variation  of  impedance  modulus  with  frequency  for  an  arbitrarily 
selected  two-degree-of-freedom  model.  Fcur  different  values  are  used  for  the  lower 
spring  and  damper.  It  is  remarkable  that,  when  the  lower  system  is  soft  such  as  in 
the  cases  shown  where  =  10  or  20,  the  impedance  curve  appears  to  be  that  of  a  single- 
degree-of-freedom  system  and  gives  no  indication  of  the  presence  of  the  second  peak  at¬ 
tributable  to  the  second-degree-of-freedom.  It  is  also  notable  that,  at  high  frequencies 

|z|  2Kj  as  £2  « 

This  result,  and  those  indicated  in  figure  1,  can  be  generalized  for  any  number  of  degrees 
of  freedom.  It  implies: 

(a)  that  the  lowest  mass  tends  to  become  stationary  at  high  frequency, 
so  that  all  the  work  is  absorbed  in  the  lowest  damper,  and 

(b)  that  the  lowest  damper  therefore  masks  the  effect  of  high  fre¬ 
quency  resonances  further  up  in  the  system. 

In  a  similar  way  it  can  be  shown  that  when  the  model  is  terminated  by  a  mass  (mn)  at  the 
driving  point  (as  in  figure  lc) 

|Z|-*-mni?  as  Q  oo 

At  the  low  frequency  end  of  the  impedance  characteristic,  the  whole  system  mass  tends  to 
move  as  a  single  body  so  that 

|z|  -*-mT£  as  £2  — o 

where  mT  js  the  total  mass.  This  result  is  true,  whatever  type  of  element  terminates 
the  model  at  the  driving  point. 

Let  us  now  pretend  that  we  do  not  know  that  the  impedance  curves  in  figure  2  were  devel¬ 
oped  by  a  two-degree-of-freedom  system.  Lacking  this  knowledge,  it  would  be  reason¬ 
able  to  assume  that  the  system  generating  the  curves  of  figure  2  could  be  approximated  by 
a  single-degree-of- freedom,  particularly  for  the  lower  values  of  Wj  .  Such  an  assumption 
would  lead  to  the  results  shown  in  table  1. 

It  is  interesting  to  note  that  this  calculation  gives  a  reasonable  figure  for  the  apparent 
stiffness  (k^)  of  the  lower  spring,  the  maximum  error  being  a  20'X  under-estimate.  The 
apparent  damping  constant  (2K^)  of  the  equivalent  single-degree-of-freedom  model  agrees 
very  closely  with  the  actual  value  for  the  lower  system  (kj),  based  on  both  Z  max  and  the 
asymptotic  value  of  Z  as  £2  oo  .  For  practical  purposes,  then,  the  apparent  damping 
ratio  (£T^)  of  the  single-degree-of-frecdom  approximation  is,  however,  considerably  less 
than  the  actual  damping  ratio  (cj)  of  the  actual  system.  Thus,  the  single-degree-of- 
freedom  assumption  gives  a  satisfactory  indication  of  the  characteristics  of  the  lower  mem 
bers  of  the  multi-degree  system. 
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The  relative  amplitudes  of  the  two  masses  are  shown  in  figure  3.  It  is  clear  that  the  Z  max 
corresponds  to  both  masses  moving  more  or  less  together  on  the  lower  spring  and  damper. 

A  second  resonance  occurs  at  a  frequency  of  about  8  Hz,  and  may  be  regarded  as  being  due 
to  the  upper  mass  moving  on  its  own  spring  and  damper,  rather  like  a  vibration  absorber, 
so  that  the  lower  mass  is  relatively  quiescent.  It  will  now  be  shown  that  this  second  reso¬ 
nance  is  independent  of  the  lower  spring  and  damper  values. 


il/2ir  FREQUENCY  (Hz) 


Values  of  2ki  -  2 


TABLE  I 


Equivalence  of  One-  and  Two-Degree -of- Freedom  Systems 


i a  (rads /sec) 

10 

20 

48.7 

70 

kj  (lb/ft) 

125.0 

500.0 

29.60 

6.20 

Values  for 

*1  <-> 

.25 

.25 

.25 

.25 

2- Degree- 

2  , 

mj  (lb  sec  /ft) 

1.25 

1.25 

1.25 

1.25 

of-  Freedom 

m^,  (lb  sec2/ft) 

3.73 

3.73 

3.73 

3.73 

System 

2K1  =  2cm1tu1  (lb  sec/ft) 

6.  25 

1.25 

30.4 

43.8 

from 

l  Z|  (lb  sec/ft) 

max  '  ' 

80.0 

158.0 

340.0 

400.0 

figure  2 

Frequency  for  |Zl 

Qo  (rad/sec)  max> 

5.98 

12.58 

27.65 

36.5 

• 

|zl  /  Q  m  •  (  -  ) 

max  o  T 

•  3.58 

3.37 

3.295 

3.24 

Apparent 

♦Apparent  c‘  (  -  ) 

A 

1.48 

.155 

.16 

.163 

Equivalent 

(.592) 

(.62) 

(.64) 

(.652) 

1-Degree- 

Apparent  u  (rad/sec) 

A 

5.98 

12.58 

27.65 

36.5 

of- Freedom 

•V"!1 

2 

(.598) 

(.  629) 

(.568) 

(.521) 

System 

Apparent  k  =  u  m 

A  T 

133.0 

469.0 

28.40 

49.60 

•  w 

(1.06) 

(.  939) 

(.  960) 

(.  81)- 

Apparent  2K  =  2c  m  u 

A  A  1  A 

(lb  sec /ft) 

6.  6 

14.5 

33.0 

44.4' 

(2Ka/2Kj) 

(1. 05) 

(1.06) 

(1.085) 

(1.01) 

*  Those  values  of  are  read  from  figure  24  of  reference  7,  for  the  appropriate  values 

of  |Z|  fQ  m  . 

1  'max  o  T 


in 


AMPLITUDE  RATIOS 


Sa=  compression  of  spring  n 
=  An*  t^n  "  ^n-1  ) 

^  unstretched  length  of  spring  n 


Figure  4.  Notation  for  a  Two-Degree-of-Freedom  System. 


The  equations  of  motion  for  the  two-degree-of-freedora  system  shown  in  figure  4  may  be 
written  as  follows:  • 


h*,  i,  *  *,0  *  to,  i 
f : i  ->  K  Z, )  J 


•  • 

•  •  c 

*  J  _  \ 


1- 


-  /V,S  •  -  <K*(-Xz  ■*  £ J  +  S, 

•  •  *# 

5”  ( ■*  to3‘i7)  A  S,  4 


where 


-v  -  c r 


C 


/  ^ 


Behavior  of  the  Lower  System 

The  force  In  the  lower  spring  and  damper  must  be  equal  to  the  applied  force 


i.e. , 


or 


2.^,  S,  —  ~~ 


K  <-  (aU*  =  " 


Fc 

ZK 


S.  XXXa*"  rCJ 


_  ej,  t 


~  JL 


*IU 


Fc  Jt  +  C 


H 


For  the  case  when  Fe  =  F0  $V*  -H-t  .  the  Integral  becomes 


X 


-e  e/f  -  3.  f. 
*  2Ki  J 


a),  t 

-O'  ftr* 


= 


^  C‘*6k3trj 


•  S'  -  JS—  “  Cctji-tl 

[•  *  (ifcjd  7 


CJ2 


fe/t  f 


(3) 


(4) 


<5) 


The  term  containing  the  integration  constant  C  is  a  transient  function  of  the  initial 
conditions,  and  can  be  ignored  when  considering  steady  state  sinusoidal  excitation,  because 


i  •« 
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■  ai  V  asidwtfterii a y  mmwWi-s  iaUirf 


■»T<  l  ■<MfW|iiiWrttfilliifiri  Irtu  I  rosiiii 


J 


■  ' 


V. 


X  ..  .x 


V.  1 


the  steady  state  condition  can  be  written  as 


U>,  -  ZC.-flCiiSLt 


[(iZJl)*  *■ 

Note  that  as  _TL  — >  O 


S,  — »  Sl»JLlT  =  &- 

*  CJ.y”,  *1 

and  yl  =  y2  =  yc 

This  Is  the  static  displacement  result,  /  £/-  £  .  where  the  force  is  transmitted  by  the 

spring  and  the  whole  system  moves  as  a  solid  body,,  in  phase  with  the  applied  force. 

When  — *  ©o 

-  — >  Q 

*  c,  . 

Fq  .  Mitt  Jt.tr  -^4  0 

a£/w 

Thus  the  motion  lags  the  applied  force  by  Tr/3.  ,  and  all  the  force  is  absorbed  in  the 

damper. 

The  maximum  value  of  S,  is  obtained  by  re-writing  equation  (6)  as 


s,  =  Fo  Sit>  (M-t  *■&,) 
OJ,*  y*},  / 1  *- 


where  (pi  = 


The  maximum  deflection  amplitude  occurs  in  the  limit  case  ( -0.-+0  )  when  all  the  force 
is  absorbed  by  the  spring  but  the  motion  amplitudes  are  so  large  when  _n.-*  ©  that  the 
displacement  amplitude  ratios  and  y,/ye  tend  to  unity,  masking  the  fact  that  the 

spring  deflection  is  maximum  at  this  point  (see  fig.  S). 


but 
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HIGHER-D.EGREE-OF- FREEDOM  SYSTEMS 

Figure  6  shows  the  Impedance  of  several  multi-degree-of-freedom  systems  with  a  low 
damping  coefficient  of  c  =  .03.  Despite  this  low  damping,  visibility  of  the  higher  modes 
is  very  poor,  the  corresponding  peaks  being  orders' of  magnitude  less  than  the  first  peak, 
and  therefore  virtually  undiscoverable  under  real  experimental  conditions. 

When  the  damping  ratio  is  increased  to  the  value  of  c  =  0.3,  more  typical  of  biodynamic 
structures,  (fig.  7)  it  is  found  that  only  the  first  two  or  three  peaks  are  visible,  further 
compounding  the  problem. 

The  corresponding  complex  impedance  plot  for  the  ten-degree-of-freedom  system  with 
c  =  0.3  is  given  in  figure  8  and  the  relative  amplitudes  of  three  of  the  ten  masses  in 
figure  9. 


IMPEDANCE  OF  A  NONLINEAR  SYSTEM 

So  far  we  have  considered- only  linear  systems.  Even  so,  only  the  first  mode  can  be  seen 
clearly  in  the  impedance  plot  when  the  damping  is  appreciable,  as  in  the  human  body. 
Phillips  (ref.  4)  found  that  force  deflection  measurements  on  a  seated  subject  obeyed  an 
approximately  cubic  deflection  law.  He,  therefore,  developed  a  single-degree-rf- freedom 
model  mounted  on  a  damped  cubic  spring.  The  cubic  spring  has  different  stiffnesses  at 
different  values  of  mean  loadings  (or  "bias")  which  conforms  with  the  observed  change  in 
impedance  characteristics  measured  in  centrifuge  experiments  (refs.  5,  7).  Jn  figure  10 
the  apparent  impedance  of  the  Phillips  model  is  presented  with  "coulomb, "  linear  and  hy¬ 
draulic  damping,  at  three  different  bias  levels. 
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MAGNITUDE  (lb  sec/ft) 


/Z/  IMPEDANCE  MAGNITUDE  (lb  scc/ft) 


for  Frequencies  From  .  01  to  100  (C  =  .  3)  f< 
Degrees  of  Freedom  Model  of  Figure  7. 


ffi 


Z  \ 

Impedance 

Magnitude 


t 


10  degrees 
of  freedom 
damping  faqtor 
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It  is  noteworthy  that  coulomb  damping  can  result  in  an  apparent  impedance  an  order  of 
magnitude  greater  than  that  of  the  same  system  with  linear  damping  and  that  a  change  in 
loading  bias  of  the  cubic  spring  can  change  the  impedance  very  considerably.  This  figure 
clearly  Indicates  that,  unless  the  type  of  damping  is  known,  i.  e. ,  viscous,  coulomb,  etc. , 
and  the  spring  deflection  characteristics,  impedance  data  tells  us  almost  nothing  about  the 
nature  of  the  system. 


THE  DEFECTS  OF  EXISTING  IMPEDANCE  DATA  FOR  HUMAN  SUBJECTS 

Impedance  measurements  of  live  human  subjects  have  been  so  far  obtained  in  two  ways  - 
analysis  of  abrupt  deceleration,  and  "steady  state"  sinusoidal  excitation.  The  most  recent 
"steady  state"  measurements,  and  in  many  ways  the  most  interesting,  are  due  to  Vogt, 
Coermann  and  Fust  (ref.  5),  (seated  position)  and  Vykukal  (ref.  6),  for  a  semisupine  position. 
In  both  cases  a  centrifuge  was  used  to  bias  the  steady  state  acceleration  and  both  investiga¬ 
tors  noted  changes  in  impedance  with  increasing  acceleration.  Some  of  this  data  is  repro¬ 
duced  in  figures  11  and  12  for  1,  2,  and  3  g  bias  conditions. 

In  both  investigations  the  experimental  scatter  appears  to  be  quite  large,  even  though  the 
experimental  points  are  (most  unfortunately)  not  given  in  the  papers,  and  can  only  be  in¬ 
ferred  from  the  plots.  This  scatter  is  not  too  surprising  when  it  is  noted  that  in  an  earlier 
study  Coermann  (ref.  7)  obtained  for  a  1  g  loading: 

|Z|max  =  6. 65  x  10®  dyne/cm/sec,  subject  sitting  erect  =  455. 5  lb  sec/ft 

A 

|Z|jnax  =  5.08  x  10  dyne/cm/sec,  subject  sitting  relaxed  =  347.5  lb  sec/ft 

At  the  low  force  levels  associated  with  shake  table  measurements,  muscle  tension  evidently 
plays  a  major  role  in  the  dynamics  of  the  body.  For  large  loadings,  where  the  failure  mode 
is  generally  vertebral  fracture,  one  would  not  expect  muscles  to  play  so  significant  a  role. 

The  most  recent  impedance  measurement  (ref.  5)  under  l  g  conditions  ("slightly  erect" 
posture)  resulted  in  an  impedance  peak  of  only  267. 5  lb  sec/ft,  possibly  indicating  an  even 
larger  inherent  scatter  than  the  two  results  mentioned  above.  A  spread  of  nearly  2:1  in 
|Z|  max  is  not  very  encouraging,  especially  when  dynamic  model  analysis  of  ejection  injur¬ 
ies  (for  example)  indicates  that  the  large  deflection  dynamic  characteristics  of  the  human 
body  are  rather  surprisingly  consistent. 

In  this  connection,  it  is  suggested  that  the  Vogt,  Coermann  and  Fust  (ref.  5)  data  is  more 
meaningfully  presented  as  shown  in  figure  12,  rather  than  by  connecting  up  the  data  points 
with  straight  lines,  and  inferring  that  the  resulting  peaks  and  valleys  have  some  physical 
meaning. 

In  a  later  section  of  this  paper,  it  will  be  suggested  that  the  first  and  largest  peak  in  the 
impedance  modulus  for  a  seated  human  is  due  to  the  whole  body  resonating  on  the  buttocks. 
This  statement  is  only  approximately  correct,  of  course,  because  other  resiliences  in  the 


23 


/z  /  IMPEDANCE  MAGNITUDE  (lb  sec/ft) 


n/2  if  FREQUENCY  (Hz) 


Figure  11.  Average  Experimentally  Measured  Impedance 
of  Seated  Hur.ipi;  Subject  Under  1  g  Conditions 
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/z  /IMPEDANCE  MAGNITUDE  (lb  sec/ft) 
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body  arc  also  deforming,  but  to  a  lesser  extent.  Since  the  variability  ir.  the  size  and  muscle 
tone  of  the  buttocks  is  greater  than  for  most  other  portions  of  the  (male)  human  body,  it  is 
not  surprising  that  impedance  measurements  dominated  by  the  buttock  mode  should  be  very 
variable,  from  subject  to  subject. 

Tlie  available  data  on  |z|max  is  given  in  table  2,  together  with  the  characteristics  of  approxi 
mately  equivalent  single-degree-of-freedom  lumped  parameter  systems  which  would  give 
the  same  value  of|'/|max  at  the  same  frequency.  The  two  values  of  1 2^  max  for  the  "sitting 
erect"  case  are  plotted  in  figure  13  and  compared  with  the  results  given  by  the  simple  two- 
degree-of-freedom  model  shown  in  figure  2,  using  the  ratios  of  the  undamped  natural  fre¬ 
quencies  calculated  in  table  2  (rather  than  their  absolute  values).  The  trend  is  seen  to  be 
in  the  right  direction,  at  least,  lending  support  to  the  hypothesis  that  Coermann's  (ref.  7) 
buttocks  in  the  1963  experiments  were  significantly  stiffer  than  the  average  of  the  ten  sub¬ 
jects  used  in  the  1968  report  (ref.  5). 

An  attempt  to  reproduce  these  measured  impedance  data  by  using  a  four-degree-of-freedom 
lumped  parameter  model  is  described  in  the  next  section. 


Experimental  Values 


TABLE  n 


Equivalent  Single- Degree-of- Freedom  Model 
for  "Steady  State"  Impedance  Measurements  at  1  g 


Sitting  Relax¬ 
ed  (ref.  7) 

Sitting  Erect 
(ref.  7) 

Sitting  ’’Slightly 
Erect’’  (ref.  5) 

lzlmax  <dWcm/sec) 

5. 08  x  106 

6. 65  x  10® 

3. 9  x  106 

(  ft  max  <lb  8ec/ft) 

347.5 

455.5 

267.5 

1  Frequency  for  |z|  max 

Qo/ in  (Hz) 

5.3 

6.3 

5.0 

/  Frequency  for  |z|  max 

Q0  (rads/sec) 

33.30 

39.6 

31.42 

I  Subject  Weight  (lb) 

185.0 

185.0 

176.0 

|  Effective  Weight*  mg 

(lb) 

161.0 

161.0 

153. 0 

Effective  Mass  ♦  m  (slugs) 

5.01 

5.01 

4.76 

/  Mmax/0°m 

(-  ) 

2.083 

2.30 

1.7884 

1  c  *♦ 

(-) 

0.27 

0.24 

0.33 

1  Oo/u  ** 

(-  ) 

1.018 

1.012 

1.037 

) 

(Hz) 

5.21 

6.22 

4.82 

1  u 

(rads/sec) 

32.8 

39.1 

30.3 

\  *  Qo/Q  res 

<  -  ) 

1.105 

1.078 

1.175 

O  /2tt 
res 

(Hz) 

4.8 

5.84 

4.25 

Where  Q  reg  Is  the  resonant  frequency  of  the  system  at  which  the  amplitude 
magnification  is  maximum. 


*  Effective  weight  is  taken  as  87%  of  the  total  weight  on 

the  assumption  that  the  lower  legs  and  feet  are  not  excited. 

**  From  Payne,  reference  8. 
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SECTION  HI 


A  LINEAR  FOUR-DEGREE-OF- FREEDOM 
SPINAL  MODEL  FOR  THE  SEATED  MAN 


The  linear  single-degree-of-freedom  "spinal"  model  currently  used  to  calculate  Dynamic 
Response  Index  (DRI)  is,  of  course,  only  the  starting  point  for  the  development  of  a  more 
sophisticated  model.  The  first  step  is  to  add  a  "pelvic  mass"  at  the  bottom  of  the  "spine 
spring"  and  a  buttock  resiliency  underneath  the  pelvic  mass.  The  pelvic  mass  and  buttock 
spring  represent  the  "buttock  mode, "  and  for  the  reasons  discussed  earlier,  it  seems  likely 
that  this  mode  has  the  greatest  influence  on  impedance  measurements  of  seated  subjects. 

The  next  logical  addition  to  the  model  is  a  degree  of  freedom  representing  the  viscera, 
since  many  studies  have  shown  that  this  has  a  marked  resonance  between  3. 0  and  3.5  Hz 
(see,  for  example,  ref.  1  and  the  quoted  exerpt  on  p.45)..  As  shown  in  figure  14,  this  is 
most  simple  shown  as  being  supported  from  the  upper  torso  mass,  although  a  case  could 
obviously  be  made  for  an  additional  spring-damper  connection  to  the  pelvic  mass,  shown 
dotted  in  the  figure.  In  view  of  the  relatively  greater  significance  of  the  links  closest  to 
the  shake  table,  it  would  appear  that  the  connection  5  will  have  a  far  greater  influence  on 
the  impedance  than  3.  However,  this  concept  (5)  has  not  been  pursued  further. 

Finally,  since  shake  table  measurements  have  shown  marked  head  resonance  with  respect 
to  the  upper  torso  at  around  30  Hz,  it  is  desirable  to  add  the  degree  of  freedom  on  top  of 
the  upper  torso. 

These  four  degrees  of  freedom  give  the  model  shown  in  figure  14.  In  the  next  section  we 
first  derive  the  linear  equations  for  this  model,  and  then  a  more  general  relationship  which 
permits  the  inclusion  of  nonlinear  terms.  The  coefficients  of  these  equations  are  then 
determined  by  an  examination  of  experimental  data,  using  various  techniques. 


s 
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EQUATIONS  FOR  THE  FOUR- DEGREE-OF- FREEDOM  MODEL 


Viscera 


Figure  14.  The  Four-Degree-of- Freedom'  Model  of  a  Seated  Man 
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X*,  =  unstretched  length  of  spring  n 

spring  deflection  (+  ^  In  compression) 

Then  =  (**»  i,  +  k>  %>)  T  (Mi  *1  *  ^  H) 

(2^i  S3  +  63X3)  ■+  (2K. +  k%&*)  —  (2^  £<,  4  k*  &+) 

(2X*%U+ kM 
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The  equations  (13)  become  ^  . 

•jt  .  fas,  ->^?0  *  "*  ^  ) 

p  -  fa  i ,  *>;  so  i  4fc)fa  J  -* '  4®1*' K  r*x) 

£  =  /24^**VtO- 

y%  =  -7 2<>  -+•  SO 
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•■••?'  ■ 
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(14) 


Making  these  substitutions,  the  equations  of  motion  (14)  become 

£  =,  Cs,  v  *■,  O  -  4(W 2'« 

g  ,  /$,  S,  ♦-;!.)  Z'*S*  y  S,y 

—  4? K  ■* n*  sv )  *  s. 

^  =  ~('i9  -*  s,  v  5, 

Note  from  (17)  and  (18)  that 

(Zu4t<*i+  *('£%' 4**^  *  ° 


(15) 

ae) 


a7> 

as) 


THE  NONLINEAR  EQUATIONS 

If  the  damping  and  spring  forces  are  nonlinear,  we  replace  the  terms, 

.  i 

,  >  {  *2c^  £  +  fj*  £^  )  by  ,  for  simplicity,  so  that  the  equa¬ 


tions  become 

•  0 

* 


■  *;  *f,  - 

■  %  +'i>  *  f>  *  4t\)h  -  W%)r„ 


a9) 


(20) 


Equations  (15)  through  (22)  show  how  the  response  of  one  degree  of  freedom  is  modified  by 
the  other  degrees  of  freedom.  For  example,  suppose  m  and  m4  are  zero  so  that  d  (2/3) 
and  <$  (4/3)  are  both  equal  to  zero  and  the  equation  (20)  becomes: 

yc  =  +  «3  +  f3 

or  for  a  linear  system  [from  equation  (17)1 

‘4  =  <*3  +  2c3*3  +  w3  S3>  +  h 
This  is  the  "cushion  theory"  equation  of  reference  6. 


THE  BUTTOCK  MODE  IN  THE  LINEAR  MODEL 

It  has  already  been  established  that  the  closer  an  element  of  a  multi-mass  model  is  to  the 
driving  point  the  more  effect  it  has  on  the  impedance  characteristics.  This  is  illustrated 
in  figures  6  and  7.  Thus,  it  is  clear  that  probably  the  most  important  components  of  the 
seated  human  body  in  influencing  the  impedance  are  the  buttocks.  It  is  also  rather  obvious, 
upon  consideration  of  this  subject,  that  the  spring  rates,  and  possibly  the  damping  rates  of 
the  buttocks  may  be  highly  nonlinear.  Under  high  g  loads,  for  example,  the  buttocks  will 
"bottom  out"  and  cease  to  influence  the  behayior  of  the  body.  The  very  markedly  different 
effects  of  various  types  of  nonlinearity  are  shown  in  figure  10.  In  order,  however,  to 
operate  the  linear  four-degree-of-freedom  model  of  the  seated  human  body,  it  is  necessary 
to  assign  linear  values  to  the  buttock  spring  stiffness  and  damping  irate  terms.  The  values 
used  were  as  follows  (the  subscript  1  refers  to  the  buttock/pelvic  mass  system): 


Pelvic  mass 

“l  = 

64/32.2 

2 

(lb-sec  /l 

Undamped  natural  frequency 

"l  = 

39.55  to  198 

(rad/sec) 

Damping  ratio 

51  -  = 

.5  to  .05 

The  subject  of  a  more  realistic  nonlinear  representation  of  the  buttocks  mode  will  be  • 
treated  in  a  later  chapter  (see  p.  88). 


THE  SPINAL  MODE  AND  ITS  APPROXIMATE  LINEAR  MODEL 

A  method  of  computing  the  stiffness  of  the  human  spine  has  been  given  by  Stech  and  Payrie 

(ref.  9).  This  method  is  based  upon  Yorra's  (ref.  11)  stiffness  measurement  for  a  single 

L4  vertebra  (fig.  15)  extrapolated  for  all  vertebrae  between  L5  and  T1  inclusive  by  means 

of  Ruff's  (ref.  2)  breaking  load  measurements  shown  in  table  3.  In  the  following  analysis 

the  subscript  3  refers  to  the  total  spinal  loads  and  deflections  and  to  the  upper  thorax  mass, 

and  the  subscript  m  refers  to  those  parameters  that  are  particular  to  one  vertebra  (m  =  L4, 

etc. ).  Subscript  1  and  2  will  be  used  to  indicate  the  deflection  of  a  particular  vertebra  from 

deflection  S  to  S 

m,  m„ 
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Calculation  of  Vertebral  Stiffness 


56*  89.6  2161*  24.1  99.6  Ruff  (ref.  2) 

**  Represents  value 

58*  .  92.8  2168  *  23.4  100. 0  assumed  by  Stech 

and  Payne  (ref.  10) 

60*  96.0  2366* *  24. 6  109.1 


The  errors  inherent  in  the  method  are  such  that  it  ts  pointless  to  strive  after  extreme 
accuracy.  Thus,  a  low  order  polynomial  fit  to  Yorra's  force-deflection  curve,  as  shown 
in  figure  15,  is  probably  satisfactory  for  deflections  less  than  one-tenth  of  an  inch.  This 
fit  given  by  an  expression  of  the  following  form  is  assumed  to  be  similar  for  each  of  the 
individual  vertebrae  with  appropriately  scaled  coefficients: 


(23) 


Where  ®m  is  the  deflection  of  the  m  vertebra.  In  the  case  of  the  L4  vertebra  shown  in 
figure  15,  the  coeffients  used  were  as  follows: 


when  24  is  in  inches 


when  is  in  feet 


=  3.88  x  10 


=  4.66  x  10 


r 


=  2. 032  x  10 


=  1.22  x  10 


6 


f 

r 


=  2.925  x  10 


=  2.11  x  10 
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The  local  stiffness  is 

Jfsj 


*  *  -r*  - 


(24) 


An  "average  stiffness"  of  a  nonuniform  spring  may  be  defined  by  equating  the  work  required 
to  compress  a  spring  with  the  characteristics  given  in  equation  (23)  (from  deflection  £», 
to  deflection  in  the  case  of  the  nth  vertebra),  with  that  required  to  compress  an  equival¬ 
ent  uniform  spring  of  stiffness  kmi2-  The  work  done  W1>  2  to  compress  each  of  these  • 
springs  is  given  by: 


Thus,  for  the  spring  represented  by  equation  (23): 

'•  -  f(«S~ys*-r£)js„ 

and  for  the  equivalent  uniform  spring  (defined  by  the  equation  F$t-  J*.  ): 

=■  ~  ) 

Equating  these  two  values  for  fe^and  making  the  simplifying  assumption  that  o  which 
can  be  done  with  little  loss  of  accuracy  then  the  average  stiffness  is  given  by: 


(25) 
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Following  the  Stech  and  Payne  procedure,  the  stiffness  between  each  vertebra  and  the  base 
of  1.5  at  1  g  lading  (table  4)  can  now  be  computed,  with  the  results  illustrated  in  figure  17. 
Stiffness  multipliers  for  different  g  values  (from  figure  15)  and  age  (from  Stech  and  Payne) 
are  given  in  figures  18  and  19. 


In  the  original  Stech  and  Payne  analysis,  the  spine  was  taken  to  extend  from  T1  to  L5  inclu¬ 
sive.  A  more  detailed  examination  of  the  human  anatomy  has  suggested  that  termination  a 
the  top  of  T4  would  be  more  logical  as  little  of  the  upper  thorax  is  supported  by  vertebrae 
above  this  point.  Thus,  from  figure  17  the  total  stiffnesses  from  L5  to  T4  are,  for  a  male 
subject  age  57.5  years 

LOCAL  1  g  stiffness  =  4.97  x  103  Ib/ft 

AVERAGE  1  g  stiffness  =  3.15  x  103  lb/ft 

Corrected  to  age  27.9  years  by  means'  of  a  multiplier  of  1.52  read  from  figure  19,  these 
give 

LOCAL  1  g  stiffness  =  7.56  x  103  lb/ft 

,  (26) 
AVERAGE  1  g  stiffness  =  4.79  x  10  Ib/ft 


For  small  amplitude  excitation,  assuming  an  upper  thoracic  mass  m3  =  1.49  lb  sec2 /ft, 
the  local  1  g  stiffness  gives 


w3  = 


/  7.51 


56  x  10" 
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=  71.1  rads /sec  (1  g) 


For  smali  amplitude  excitation  about  2  g  and  3  g  we  have,  from  figure  18 

LOCAL  2  g  stiffness  =  7.56  x  1.255  x  103  =  9.50  x  103  lb/ft 

3  g  stiffness  =  7.56  x  1.43  x  10  =  10.80  x  10  lb/ft 


and  w3  =  71.1  x  /l. 255  =  79.6  rads/sec  (2  g) 

=  71.1  x  yi.43  '  =  85.0  rads/sec  (3  g) 


(27) 

(28) 


These  stiffnesses  and  frequencies  are  substantially  greater  than  the  figures  previously  de¬ 
duced  by  Payne.  At  10  g  the  average  stiffness  is  given  by 


Average 
**>3  = 


V8.9 
— 


10  g  stiffness  =  4.79  x  1.855  x  103  =  8. 9  x  103  lb/ft 


xJO 
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-  =  77.4  rads/sec  (10  g) 


(28) 


Corrected  to  the  single-degree-of-freedom  spinal  model,  where  the  head  (12  lb)  and  upper 
torso  (48  lb)  masses  are  lumped  together,  we  obtain 


w3  =  77 


4  x 

•  y  6o.o 


=  69.3  rads/sec  (10  g) 


(30) 


=  11.05  Hz  (10  g) 


TABLE  IV 


1  g  Stiffness  Calculations  flb/ft  units) 
(for  57.5  Year  Old  Male  Subject) 


Stiffness  of  Each  Vertebra 

Total  from  L5  to  Each  Vertebr^ 

% 

Vertebra 

Local  Stiffness 

Average  Stiffness 

Local 

Average 

.166 

T1 

18. 02  x  103 

(from  0  to  lg) 
11.4  x  103 

3.016  x  103 

1. 910  x  103 

.221 

T2 

24.0 

15.2 

3. 623 

2. 295 

.277 

T3 

30.1 

19. 05 

4.27 

2.703 

.332  • 

T4 

36.0 

22.8 

4.972 

3.149 

.387 

T5 

42.0 

26.6 

5.768 

3.654 

.461 

T6 

50.0 

31.7  • 

6.  69 

4.236 

.535 

T7 

58.1 

36.8 

7.72 

4.89 

.607 

T8 

65.9  • 

41.7 

8.90 

5.64 

.689 

T9 

74.8 

47.4 

10.29 

6. 52 

.753 

T10 

81.7  * 

51.8 

« 

11.93 

.7.56 

.784 

Til 

85.1 

53.9 

13.97 

8.85 

.810 

T12 

87.9 

55.7 

16.72 

10.59 

.826 

LI 

89.6 

56.8 

20.  6b 

13.08 

.888 

L2 

96.4 

61.0 

26.  83 

16.  99 

.996 

L3 

108.2 

68.5 

‘37.18 

23.55' 

1.0 

L4 

108.5 

68.7 

56.  64 

35.88 

1.091 

L5 

118.5 

75.1 

118.5 

75.1 
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m  VERTEBRA 


Figure  17.  Spinal  Stiffness  at  1  g  for  57. 5  Year  Old  Male 


EQUIVALENT  STEADY  STATE  ACCELERATION  (g) 


Figure  18.  Average  and  Local  Stiffness  Between  L5  and  T4 

Inclusive  (Age  57.5  Years)  as  a  Function  of  Loading 


AGE  (years) 


Figure  19. 


Correction  Factor  for  the  Effect  of  Age  on 
Spinal  Stiffness 


The  equivalent  values  using  the  T1  vertebra  instead  of  the  T4  vertebra  would  be: 
AVERAGE  10  g  stiffness  (Tl)  =  1. 91  x  103  x  1.52  x  1.  855  =  5. 39  x  103  lb/ft 

8Ild  /5i39  x  103  x  32.2  ..  ...  . 

T1  V - ~60 -  =  54  rad//sec  (10  8) 

Cj3(TD/2^  =  8.58  Hz 


This  is  very  close  to  the  value  of  8.45  Hz  deduced  for  the  original  Payne  injury  model. 

However,  it  is  now  assumed  to  be  more  realistic  to  use  the  T4  values  for  k„  and  u  . 

3  3 

In  the  linear  four-degree-of-freedom  model  described  in  the  next  chapter  of  this  report, 
the  following  ranges  of  values  were  chosen  for  a  parametric  study  of  the  impedance  of 


spinal/upper  torso  mode: 

Upper  Torso  Mass: 

“3 

= 

48/32.2  lb  sec2 /ft 

= 

1.49  lb  sec2/ft 

Spring  Rates: 

k3 

= 

2400  to  6300  lb/ft 

Undamped  Natural  Frequencies: 

=  • 

40  to  65  rads/sec 

Damping  Ratio; 

*3 

= 

.075  to  .3  (  -  ) 

Damping  Rates: 

2K 

= 

13. 2  to-52.  8  lb  sec/ft 

These  values  were  varied  in  the  ranges  given  in  an  attempt  to  match  measured  impedance 
curves  such  as  those  shown  in  figures  ?i  and  12.  The  mass  was  fairly  readily  definable 
and  so  was  not  varied.  The  spring  rate  k  was  varied  across  a  range  representing  the 

spread  of  observed  data.  The  damping  was  varied  over  a  much  wider  range  both  to  ob¬ 
serve  the  effect  of  the  change  and  because  experimental  measurements  vary  widely. 


THE  VISCERAL  MODE 


Goldman  and  Von  Gierke  (ref.  1 )  have  pointed  out  that  "one  of  the  moat  important  sub¬ 
systems  of  the  body,  which  is  excited  in  the  standing  and  sitting  position  as  well  as  in  the 
lying  position,  is  the  thorax-abdomen  system.  The  abdominal  viscera  have  a  high  mobil¬ 
ity  due  to  the  very  low  stiffness  of  the  diaphragm  and  the  air  volume  of  the  lungs  and  the 
chest  wall  behind  it.  Under  the  influence  of  both  longitudinal  and  transverse  vibration  of 
the  torso,  the  abdominal  mass  vibrates  in  and  out  of  the  thoracic  cage.  Vibrations  take 
place  in  other  than  the  (longitudinal)  direction  of  excitation;  during  the  phase  of  the  cycle 
when  the  abdominal  contents  swing  toward  the  hips,  the  abdominal  wall  is  stretched  out¬ 
ward  and  the  abdomen  appears  larger  in  volume;  at  the  same  time,  the  downward  deflec¬ 
tion  of  the  diaphragm  causes  a  decrease  of  the  chest  circumference.  At  the  other  end  of 
the  cycle  the  abdominal  wall  is  pressed  inward,  the  diaphragm  upward  and  the  chest  wall 
is  expanded.  This  periodic  displacement  of  the  abdominal  viscera  has  a  sharp  resonance 
between  3  and  3.5  Hz.  The  oscillations  of  the  abdominal  mass  are  coupled  with  the  air 
oscillations  of  the  mouth-chest  system. " 

Weis  and  Mohr  (ref.  12)  have  used  cineradiographic  analysis  to  record  the  dynamic  de¬ 
flection  of  visceral  components  during  and  after  a  positive  spinal  impact  in  the  seated 
position.  They  show  the  deflection  time  history  of  the  diaphragm  as  a  typical  example. 

U  this  is  regarded  (for  the  moment)  as  a  single-degree-of-freedom  system,  then  the 
apparent  damping  and  frequency  can  be  determined  from  the  overshoot  above  the  stat¬ 
ically  deflected  position  or  datum,  measured  from  the  maximum  downward  deflection 
position,  using  the  equation 


„  .  .  „  ..  oversnoot  above  datum 

Overshoot  Ratio,  OR  =  - — -  =  1  +  e 

datum 

S,  /r 

(31) 

where  »,  *  ,J  |-  Z*  and  subscript'  2  refers  to  the  visceral  system. 

tr ' /c  = 

(The  relationship  between  c^ 

and  d  // is  plotted  in  figure  20.) 

When  this  is  done,  the  following  results  are  obtained: 

• 

Subject  Tense 

Subject  Relaxed 

Overshoot  Ratio,  OR-  (  -  ) 

1.063) 

1.212] 

1  from  ref. 

Apparent  frequency  (Hz) 

/from  ref.  12 

1 

j 

9.-68  ) 

6.37  ] 

1 

V  L  A)  = 

1 1 

2.77 

1.554 

c2  (from  fig.  20)  (  -  ) 

0. 66 

0.44 

Undamped  natural  frequency 
*./;?»/  '2r-v“^F; )  (Hz)1 

*  12.0 

7.1 

The  apparent  stiffness  is  therefore  much  greater  than  is  indicated  by  the  shake  table  mea¬ 
surements  of  Coermann  et  al  (ref.  7),  indicating,  at  first  sight,  considerable  nonlinearity 
in  the  visceral  springs.  On  the  other  hand,  the  apparent  frequencies  are  of  the  same  order 


♦See  figure  22  of  ref.  8 
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as  the  spinal  mode,  so  that  the  Weis  and  Mohr  measurements  of  the  diaphragm  motion 
may  be  showing  the  "whole  body"  response  with  some  modulation  due  to  visceral  move¬ 
ments  with  respect  to  the  rest  of  the  (deforming)  body.  In  this  report  the  latter  inter¬ 
pretation  is  assumed  to  apply.  This  means  that  the  frequency  of  the  visceral  mode  must 
be  taken  from  the  small  amplitude  measurements  of  Coermann  et  al  (ref.  7),  and  the 
value  initially  selected  is 

(Jo  -  20.4  rads/sec  (31A) 

For  an  assumed  visceral  weight  of  =  15  lbs.  this  corresponds  to  a  spring  rate  given 
by: 

k2  w22m2=  20. 42  x  =  194  lb/ft  (31B) 

There  is  no  simple  and  direct  way  to  determine  what  damping  factor  to  use.  If  the  value 
used  is  too  low,  however,  the  visceral  mode  shows  up  as  a  peak  on  the  impedance  plot. 

No  such  peak  is  observed  in  the  experimental  data.  Thus,  the  visceral  damping  is  ad¬ 
justed  to  the  minimum  value  needed  to  surpress  the  visceral  peak  in  the  impedance  magni  ¬ 
tude;  values  ofc  of  .125,  .25  and  .5  were  used.  Clearly,  the  representation  of  the  viscera 
as  being  simply  suspended  from  the  upper  torso  is  a  gross  oversimplification  as  there  must 
be  a  considerable  visceral  load  acting  directly  on  the  pelvic  mass. 

Thus,  for  the  visceral  mode  the  following  values  were  adopted  for  this  study.- 


Mass: 

m2  = 

15/32.2  lbsec2/ft 

Spring  Rate; 

k2  = 

194  lb/ft 

Undamped  Natural  Frequency: 

W2  = 

20.4  rads/sec 

Damping  Ratio: 

F2  = 

.125,  .25,  .5  (  -  ) 

Damping  Rate: 

2K2  = 

2.375,  4.75,  9.5  lb  sec/ft 
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THE  NECK  MODE 


The  transmissability  experiments  of  Dieckmann  (ref.  13)  indicate  a  head/shoulder  reson¬ 
ance  of  3.5  at  30  Hz,  for  a  seated,  sinusoidally  excited  subject, 

l.e.,  «  3.5 

)max 

If  the  motion  is  sinusoidal.  The  frequency  of  30  Hz  is  so  much  greater  than  the  values  de¬ 
rived  for  the  visceral  and  spinal  modes  that  preliminary  calculations  can  be  based  upon 
the  simplified  two  mass  system  of  figure  21. 


~  a  ^  K  +KU)  j 

\ 

/ 


And  since  y^  =  y'g  -  $ ^  the  second  equation  becomes 


-  Sv  4  2r^  -t  ^  c,  5. 
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*■  *  ci['  *  @(i)l  K  *■&*[•  +  <P(fjj  -  <p(i)£e  *>•  -A-* 


where 


m  =  ^ 


Thus,  the  apparent  undamped  frequency  Is 


CJAtt  -  -//#-  $($) 


and  the  apparent  damping  coefficient  is 


-  c,  h  +  4>(fi} 


The  apparent  damping  ratio  is 


~  ~  c*  S^W 

The  response  of  equation  (34)  will  be 

S,  =  ^7  ^  B51 

*  _*/ 

•> '  .  .  /i  Cr  /  n.  v*!1  —  a  /  _/J-  \>)  * 


s  =  ib-mT 
•' ^  f»,  -  <P(‘)  r*„  u 


where  <T«, 


is  the  steady  state  deflection  of  the  neck  spring. 


This  solution  may  be  expressed  in  terms  of  (y4)max/  (y*3)max  and  since,  from  equations 
(33)  and  (35) 


D*  <P(i)] ft 


1 


f/ 


Similarly,  since  =  2f,/}  *  cjj  S*  from  equation  (32),  and,  from  equation  (35) 


and 


_  *<j)n  t-.  m 

>»7, 


fr>  ^  c„(sit*</'J+  Si* (jit +  <yjj 


W-.  •  f^fkj  A  r«/^W 


m 


no  r 
■mox 


(37) 


=  y^w 


(38) 


This  function  is  discussed  in  Appendix  I,  where  it  is  thought  in  figure  1, 2  that  the  amplitude 
ratio  of  3. 5  corresponds  to  a  damping  coefficient  ratio  of 


c^  =  0. 15 

and  to  resonant  frequency/natural  frequency  ratio  of 
.t£b*S  -  .77? 

Ci/ty 

Thus,  if  resonance  occurs  at  30  Hz  the  frequency  cJ f  is  given  by: 

u.  =  2M  =  30.6  Hz  =  192.3  rads/sec 
4  0. 979 


(39) 


(40) 


and  the  corresponding  damping  force  is: 


—2 .  c,  cJi)  ^  —  - 


32 .17 


-  to.  7 ar  ^ 


(41) 
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Taking  a  typical  head  weight  of  12. 0  lbs,  the  spring  rate  corresponding  to  (J^  is  therefore 


k4  =  %m4 


Ifyl  x  (192. 3)2  =  14,  000  lb/ft 


Thus,  for  the  head  and  neck  mode  the  following  values  were  selected: 

Mass:  m.  12/32.2 ’.b  sec2/ft 

4 

Spring  Rate;  =  14, 000  Ib/ft 

Undamped  Natural  Frequency:  o>.  =  192.3  rads/sec  (43) 

4 

Damping  Ratio:  c.  =  0. 15 

4 

Damping  .Rate:  2K.  =  10. 75  lb  sec/ft 

4 


These  quantities  were  not  varied  in  the  subsequent  analysis  as  they  seem  to.be  fairly  well 
established  (in  ref.  13)  and  the  head  and  neck  are  remote  enough  from  the  driving  force 
for  their  influence  on  the  measured  impedance  to  be  small. 


IMPEDANCE  CALCULATIONS  FOR  AN  EQUIVALENT  LINEAR  MODEL 

Having  now  an  approximate  idea  of  the  magnitudes  of  most  of  the  coefficients  for  the  four- 
degree-of-freedom  model,  it  is  instructive  to  determine  the  driving  point  impedance  of  the 
model  and  compare  the  results  with  the  existing  experimental  data.  By  systematically 
changing  the  model  coefficients,  it  is  then  possible  to  obtain  better  agreement  with  exper- 
ment,  and  to  study  the  influence  of  the  individual  elements.  This  has  been  done,  and  some 
of  the  more  interesting  results  are  presented  in  figures  22  and  51  inclusive.  The  effect  of 
the  head  and  visceral  masses  on  the  driving  point  impedance  could  only  be  minor,  so  these 
masses,  springs  and  dampers  were  not  varied  in  this  analysis,  except  for  the  visceral 
damping  discussed  earlier.  Neither  was  the  mass  distribution  varied  as  this  was  considered 
to  be  well  enough  established. 

The  fixed  values  were  the  following 


Masses: 

Pelvic 

mi 

= 

64.0/32.2  lb  sec^/ft 

Viscera 

X 

= 

15.0/32.2  lb  sec^/ft 

Upper  Torso 

Ct 

% 

= 

48.0/32.2  lb  seCg/ft 

Head 

<■ 

= 

12.0/32.2  lb  sec  /ft 

Springs; 

Viscera/Upper  Torso 

= 

20. 4  rads/sec  (seeeq.31A) 

Head/Upper  Torso 

L 

= 

192.3  rads/sec , 

< 

= 

14, 000  Ib/ft  <6eee(l-41) 

Dampers: 

Head/Upper  Torso 

ff4 

= 

.  15  (see  eq.  39) 

and  the  variable  values  were  varied  according  to  table  5. 

Figures  22  to  51  are  plotted  to  the  same  scale  as  the  experimentally  measured  values  shown 
in  figures  11  and  12.  The  closer  t'.at  the  analytical  curve  simulate  the  measured  values  of 
figures  11  and  12,  the  more  representative  the  analytical  model  may  be  assumed  to  be. 

The  curves  plotted  were  calculated  using  an  analog  simulation  routine  on  a  digital  computer. 
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TABLE  V 

Variables  U3ed  in  Linear  Four- Degree-of-  Freedom  Model 


Mode 

2.  Viscera 

3.  Spine 

1.  Buttocks 

Link 

Damper 

Spring 

Damper 

Spring 

Damp* 

Parameter 

c2 

^3 

°>1 

ll 

Units 

rad/sec 

-- 

rad/sec 

— 

Case  No.  2 

Fig  22 

.25 

59.1 

.15 

— 

— 

Case  No.  3 

Fig  23 

.125 

59.1 

.075 

39. 55 

.125 

Case  No.  4 

Fig  24 

.25 

59.1 

.15 

39.55 

.25 

Case  No.  5 

Fig  24 

.25 

59.1 

.15 

47.5 

.25 

Case  No.  6 

Fig  25 

.25 

59.1 

.15 

55.5 

.25 

Case  No.  7 

Fig  25 

.25 

59.1 

.15 

63.2 

.25 

Case  No.  8 

Fig  26 

.25 

59.1 

.15 

198.0 

.05 

Case  No.  9 

Fig  27 

.25 

50.6 

.15 

39.55 

.25 

Case  No.  10 

Fig  28 

.25 

50.6 

.15 

39.55 

.5 

Case  No.  11 

Fig  29 

.5 

59.1 

.15 

39.55 

.25 

Case  No.  12 

'  Fig  30 

.25 

59.1 

.3 

'  39.55 

.25 

Case  No.  13 

Fig  31 

.5 

59.1 

.3 

39.55 

.25 

Case  No.  14 

Fig  32 

.25 

59.1 

.15 

47.5 

.5 

Case  No.  15 

Fig  33 

.25 

50.0 

.15 

47.5 

.5 

Case  No.  16 

Fig  34 

.25 

40.0 

.15 

47.5 

.5 

Case  No.  17 

Fig  35 

.5 

59.1 

.15 

47.5 

.5 

Case  No.  18 

Fig  36 

.5 

59.1 

.15 

44.0 

.25 

Case  No.  19 

.Fig  37 

.5 

59.1 

.15 

44.0 

.5 

Case  No. 20 

Fig  38 

.5 

59.1 

.15 

47.5 

.5 

Case  No.  21 

Fig  39 

.5  • 

59.1 

.15 

51.0 

.5 

Cas^  No.  22 

Fig  40 

.5  . 

59.1 

.15 

47.5 

.25 

Case  No.  23 

Fig  41 

.5 

59.1 

.15 

47.5 

.125 

Case  No.  24 

Fig  42 

.5 

59.1 

.15 

— 

— 

Case  No. 27 

Fig  43 

.5 

65.0 

.15 

—  . 

— 

Case  No.  28 

Fig  44 

.5 

65.0 

.125 

— 

— 

Case  No.  29 

Fig  45 

.5 

65.0 

.1 

— 

— 

Case  No.  30 

Fig  46 

.5 

59.1 

.11 

47.5 

.25 

Case  No.  31 

Fig  47 

.5 

59.1 

.15 

108.0 

.25 

Case  No.  32 

Fig  48 

.5 

59.1 

.15 

80.0 

.25 

Case  No.  33 

Fig  49 

.5 

59.1 

.15 

120.  0 

.25 

Case  No.  34 

Fig  50 

.5 

59.1 

.15 

140.0 

.25 

Case  No.  35 

Fig  51 

.5 

59.1 

.15 

160.  0 

.25 

MAGNITUDE  (lb  sec/ft) 


/Z  /  IMPEDANCE  MAGNITUDE  (tb  sec/ft) 


/  IMPEDANCE  MAGNITUDE  (lb  sec/ft) 


/z  /IMPEDANCE  MAGNITUDE  (lb  sec/fg 


/z  /  IMPEDANCE  MAGNITUDE  (lb  sec/ft) 


12.0  K'4  Head 

lbs 

19^3*5^* ',5  Nock 


48.0  lbs  m3  Uppcr 

agrf  4M  w- 

«*lM]  £&  T3‘  •“ 

©rjr  pelvic 

L^°.|m‘  M..S 

39.55?  ^ 

Buttocks 


5  10  '  15  20 

n/2  7T  FREQUENCY  (Hr) 


Figure  27.  Four  Mass  Systems  Case  No.  9 


/Z  /  IMPEDANCE  MAGNITUDE  _  (lb  sec/ft) 


/  z  /  IMPEDANCE  MAGNITUDE  (lb  sec/ftt 
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m2 _  59?I>  T  °3  ”  *15 


t.O  !;n 


Pelvic 
1  Mass 


6>]  =S 

4 


'S“l  =  .5 
"™~  Buttocks 


Figure  32.  Four  Mass  Systems  Case  No.  14 
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Figure  34.  Four  Mass  Systems  Case  No.  16 
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Figure  37.  Four  Mass  Systems  Case  No.  19 
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Figure  42.  Four  Mass  Systems  Case  No.  24 
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/%  /  IMPEDANCE  MAGNITUDE  (lb  scc/ft) 


/  z  /  IMPEDANCE  AMPLITUDE  (lb  sec/ft) 


/z  /  IMPEDANCE  AMPLITUDE  (lb  sec /ft) 


80 


SJCY  (Hz) 


I 


•*****.•« 


•  v *  ;■•-  •  ..  •  - 

^  ;„•  ;,  /'  l;  J.;;i 

■  ;  ‘  -ioi^xf-:  •■  V  ■ 

,,  \. ■■>■•'  ...  ;  •  .. 

**.  J:f-'  ^•'••,*.  •»•*«  •’  ^  •,*;•-;/•»'•'•  'y*  •;  ‘."is  •'■-..•  .  £  *  .  *  ... 

.•  .  .  „••  w  ••'  V/  'l  «V'  i,*  i.  ."»  •?•  i<'-.  •'.  ..  *;•' i  •?••  :•  .  *v  ••' . 

••' ;; ■?: .-;  «^CTION  IV 


?f%  '(£$ '  i’’ 'v  t  • ; 


■i  ?,i-"  -.  '.i; :■ . '  •"*'•  j"" 

■  tf-  .lS'r ;; 


-  •  -:  DEVELOPMENT  OF  %*. 

THE  NON-LINEAR  FOUK-DEGREE-OF-FREEDOM  MODEL 


:A v;:'^ " 

s-  -  ■  c-  '■  * ! 


The  linear  model  described  in  the  preceeding  pages  had  some  evident  shortcomings,  par¬ 
ticularly  in  the  attempt  to  force  the  spinal  and  buttock  modes  Into  a  linear  representation, 
when  both  are  clearly  highly  nonlinear.  The  next  stage  In  the  development  was  to  expand 
the  computer  model  to  accept  nonlinear  representations  of  at  least  the  spinal  and  buttock 
spring  and  damping  rates.  This  analysis  Is  discussed  in  the  following  pages. 


CORRECTIONS  TO  SPINE  STIFFNESS  AND  DAMPING  BASED  ON  IMPEDANCE  DATA 

Steady  state  mechanical  impedance  measurements  of  seated  subjects  have  been  made  under 
1  g  conditions  by  several  investigators.  As  indicated  in  figure  11,  the  results  are  very 


variable.  A-  . 

-■  •  ■  ,-V 

Vogt,  Coermann  and  Fust  (ref.  5)  also  measured  Impedance  in  a  centrifuge,  with  steady 
state  acceleration  biases  of  2  and  3  g.  Their  results  are  plotted  as  experimental  points  in 
figure  12,  because  there  seems  to  be  no  justification  for  drawing  lines  from  point-to-point 
when  the  1  g  data  has  already  demonstrated  rather  extreme  variability,  and  therefore, 
presumably,  considerable  experimental  scatter. 

The  most  reliable  1  g  data  is  presumably  the  most  recent  (shown  as  curve  3  in  figure  11). 
Figure  46  shows  a  linear  dynamic  model  which  gives  good  agreement  with  measured  data, 
up  to  a  frequency  of  about  8  Hz.  In  this  case  the  location  and  magnitude  of  the  first  peak  ts 
largely  dependent  upon  the  stiffness  and  damping  of  the  buttock  mode.  At  the  present  time  the 
poor  agreement  above  8  Hz  is  not  felt  to  be  significant,  because,  as  already  shown,  this  is 
characteristic  of  nonlinear  systems,  and  the  buttock  spring  is  known  to  be  very  nonlinear. 

It  is  also  known  that  the  buttocks  tend  to  vary5  from  subject  to  subject,  more  than  any 
other  part  of  the  human  anatomy.  Thus,  the  primary  dependence  of  the  first  peak  on  the 
buttock  mode  elegantly  explains  the  wide  variability  In  the  1  g  measurements. 

It  might  be  thought  that  the  larger  Impedance  peak,  at  a  higher  frequency  (about  8. 5  Hz)  • 
which  is  measured  under  2  and  3  g  bias,  is  caused  by  the  stiffening  of  the  buttocks  under 
the  higher  loading.  As  indicated  in  figure  52,  this  turns  out  not  to  be  so.  Increasing 
buttock  stiffness  beyond  a  value  corresponding  to  about  O  j  a  60  rads/sec  causes  the  spinal 
resonance  mode  to  increase  in  amplitude,  and  the  first  peak  becomes  asymptotic  to  about 
8. 8  Hz  as  Fo»  large  values  of  6^  the  motion  is  almost  entirely  In  the  spinal 

mode,  and  the  buttock  mode  can  be  considered  to  be  "shorted  out". 

With  "shorted-out"  buttocks  the  model  shown  in  figure  45  gives  good  agreement  with  Im¬ 
pedance  measurements,  below  about  11  Hz.  Note  that  In  figures  45  and  46 

o>3  *  59. 1  rads/sec  for  1  g,  (k„  *  (59. 1)2  x  =  5200  lb/ft)  (44) 


o> j  »  59. 1  rads/sec  for  1  g, 


<cq  « \2  v  48«  0 

<59a)  50 


WimM 


vs 

■ 


OJ3  m  65.0  rads/sec  for  2-3  g,  (kg  -  (65. 0)2  x~~-  *  6300  lb/ft) 
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BUTTOCK  FREQUENCY  «  (RADS/SEC) 


Figure  52.  Influence  of  Buttock  Stiffness  on  Location  of  the  First 
and  Second  Impedance  Peaks 


The  1  g  stiffness  given  above  is  83%  of  the  spinal  model,  and  the  2-3  g  local  stiffness  is 
79%  of  the  spinal  model,  corrected  to  age  27. 9  years.  An  average  stiffness  reduction 
to  81%  of  the  age  27. 9  years  value  is  seen  from  figure  19  to  correspond  to  a  test  sub¬ 
ject  age  of  about  44  years,  which  is  not  impossible,  but  is  unlikely.  (Unfortunately,  sub¬ 
ject  age  is  not  usually  documented  in  such  test  reports. )  On  balance,  it  seems  best  to 
provisionally  assume  that  the  live  human  spinal  stiffness  is  80%  of  the  value  determined 
from  cadaver  data. 


In  order  to  estimate  the  total  spine  force  deflection  characteristics,  using  the  data  shown 
in  figure  16  for  the  L4  vertebra,  the  deflection  scale  was  assumed  to  be  inversely  pro¬ 
portional  to  the  stiffness.  Thus,  a  deflection  of  .  1  inches  of  the  L4  vertebra  was  assumed 
to  correspond  to 

—  x  10.8.- 5  =  .1196  feet 
12  7.56 

of  the  whole  spine  between  the  L5  and  T4  vertebrae.  ( 108. 5  x  103  lb/ft  is  the  local  1  g 
stiffness  for  the  L4  vertebra  from  figure.  16,  7. 56  x  103  is  the  total  local  1  g  stiffness 
for  the  spine  of  a  27. 9  year  old  man,  from  equation  (26).]  Reducing  the  force  scale  by  a 
factor  of  .  80,  as  suggested  above,  allows  the  force  deflection  curve  for  the  whole  spine 
of  a  live  27. 9  year  old  man  to  be  given  by  making  the  appropriate  corrections  to  the  co¬ 
efficients  of  equation  (23). 

\  -  3  88  «  103  (. 8  X  )  S3  *2.032  x  105  (.«>  .TKf  »2I  S32- 

- .1.22  X  106  ,.8X  )  S33 

=  2.59  x  103  $  +  1. 138  x  105  S*  -  5. 7  x  105 C  3  (for  £  <.  12)  (46) 

3  3  v  3  3 

•  • 

In  practice  this  curve  is  found  to  produce  sensible  values  of  for  S3  4  .  12  ft.  For 
S^>  .12  ft  the  tangent  to  the  curve  at  that  point  was  used  in  numerical  calculations: 

F$3  =  331.2  +  527.8  £  3  (for  Sg^..l2)  (47) 

This  revised  curve  of  force  Fg.  versus  deflection  for  the  whole  spine  is  shown  in  figure  53. 

So  far  as  spinal  damping  is  concerned,  a  value  of  c  =  0. 1  (associated  with^g  =  65  rads/sec) 
seems  to  give  best  agreement  with  the  centrifuge  impedance  data.  Thus,  the  damping  fore.e 
is  given  by: 

FD3  =  2x0.  lx  65.  Ox  S3  =  19  483  lb.  (48) 

In  the  absence  of  any  other  data,  this  damping  is  assumed  to  be  linear  and  thus  applies 
to  the  spinal  model  for  all  deflections  and  deflection  rates. 
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BUTTOCK  STIFFNESS  f 

f  '  _  . 

As  shown  by  Phillips,  *  as  indeed,  would  be  intuitively  supposed,  buttock  stiffness  is 
nonlinear,  at  least  up  to  1  g  loading.  The  Hiillips  data,  plotted  in  figure  54  is  very 
scattered,  and  does  not  go  above  the  lg  loading  case.  One  might  suspect  that  creep 
(primarily  displacement  of  blood)  and  muscle  tension  variation  would  both  contribute  to 
this  scatter. 

An  average  value  of  Phillips*  data  is 
Fgj  =  77.0  $,  (inches)2,44 

-  3.3  S^-^xlO4  <«> 

when  S,  is  in  feet. 

The  local  slope  is 

=8.05  1,44  x  104 

d 

a  eq 

=  172.0  Fgx  ftb/ft)  (50) 


If  all  the  subject  above  the  buttock  cushion,  including  thighs,  were  rigid,  the  small  pertu- 
bation  natural  frequency  for  a  172-lb  subject,  with  150  lbs  supported  on  the  buttocks  and 
68-lb  pelvic  mass  would  be  given  by: 


/m.) 

1  32.2  . 


M0  =  26.63  (radians/sec)  (51)  • 

where  N  =  steady  state  acceleration '  (_) 

g 

Thenatural  frequency  of  the  buttock  spring  at  the  corresponding  steady  state  acceleration, 
when  associated  only  with  the  pelvic  mass  is  given  by: 


=  1-485  ^ 


(52) 


♦Unpublished  trial  report  by  Norman  J.  Phillips  on  Contract  No.  AF  33(657)-1894,  Aerospace 
Medical  Research  Laboratory,  Wright-Patterson  Air  Force  Base  (Nov  1967)  entitled  "Research 
on  Human  Responses  to  Complex  Vibrations  and  Design  Principles  for  Body  Support,  Restraint 
and  Vibration  Isolation  Systems" 


N 

1 

2 

3 

N.295  = 

1.0 

1.227 

1.383 

0^ 

II 

26.63 

32.65 

36.8 

rads/sec 

= 

4.24 

5.2 

5.86 

Hz 

<*1 

3^.6 

48.5 

54.7 

rads/sec 

In  order  to  obtain  agreement  with  impedance  measurements,  it  is  necessary  to  assume 
ui,  =  47.5  rads/sec(l  g)  and  cj,  ^  160  rads/sec  (2-3  g).  (Compare  figures  42  and  51, 
for  example.)  The  1  gvalue  of  47.5  rads/sec  is  120%  of  the  value  indicated  by  Phillips' 
data,  a  difference  which  is  not  at  all  surprising  when  it  is  recalled  that  buttock  size  and 
muscle  tone  are  very  variable.  Because  of  the  poorness  of  this  data,  therefore,  it  is  con¬ 
venient  to  adopt  Phillips’  assumption  that  the  buttock  spring  is  a  cubic  up  to  1  g. 


Rough  measurements  made  by  loading  up  a  subject  seated  on  a  hard  chair,  indicate  that  the 
buttocks  "bottom  out"  at  1.5  g,  so  that  above  this  force  value  the  buttock  frequency  is 
arbitrarily  assumed  to  be  u),  -  160  rads/sec  corresponding  to  a  "stiffness"  kj  given  by: 


k.  =  1602  x  64 
1  32.2 

=  50,100  lb 
it 


Below  1.5  g,  if 


Jfi. 


kC 


(53)  • 


At  1  g,  assuming  a  subject  with  139  lb  supported  on  the  buttocks  and  a  64.0  lb  pelvic  mass: 


t  df%,  jL  _ 

cH, 


*A 


* 
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For  the  1  g  case  where  =  1391 0  lb 

u)t  -  I*?1  ~  4>*~ 


(54) 


(55) 


and  the  general  expression  for  F^  for  less  than  a  1.5  g  load,  which  corresponds  to  a 
deflection  of  0. 1066  ft  is: 


FSi  =  1.72  x10s  £*  (56) 

Thus  the  buttock  spring  is  assumed  to  be  as  shown  in  figure  55. 

The  buttock  damping  is  assumed  to  be  constant,  so  that 

2K  =  2cm  =  2cWm  =  2  x0.25  x  47.5  x  ffLQ  =  47.2  lb  sec/ft  (57) 

1  32.2 


and  the  buttock  damping  force  is  therefore  given  by: 

Fj>|  =47.2S,  (lb)  '  (58) 


FINAL  MAN  MODEL 

The  final  form  mass  analog  of  a  seated  man  used  in  subsequent  computations  is  shown  in 
figure  5.6.  The  total  weight  of  the  man  is. assumed  to  be  160  lb,  distributed  as  follows: 


Head 

12  lbs 

Upper  Torso 

48  lbs 

Viscera 

15  lbs 

Lower  Pelvic  Mass 

64  lbs 

. 

739  lbs 

Legs  and  Feet 

21  lbs 

160  lbs 

(supported  by  buttocks) 
(assumed  to  move  with  seat) 


The  values  used  for  the  spring  and  damping  forces  for  the  four  modes  are  shown  in 
figure  56.  The  spring  and  damping  constants  art  shown  in  table  6. 
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°S1  DEFLECTION  (ft) 


Figure  55.  Buttock  Spring  Force/Deflection  Curves 

Assumed  for  a  27.  9  Year  Old  Live  Male  Subiect 


FD  =  2x0. 15x  192.3x12.0  8, 
4  32.17  4 


8  Positive  in  Compression 


Fs  =  192.3  .x  12. Ox  8a 
4  32. 17 


W3  =  48. 0  lbs. 


=  2x0.5x20.4x15.0x8,, 
!  32.17 


FSo  " 


’2  <20^15. 

f  32.17 

W«  =  15.0  lb.  „ 

<  ss 


J  % 

19.4  8 


5.3xl04xl7.2  for  0  ^$0. 1066 
0  for  8  j  <  0 

208. 3522+5. 09xl04  (8,-  0.1066) 
for  8^  0. 1066  |  W5  =  53"  5" 'll: 


17  n 

r  < 

ss 

> 

> 

S3 

> 

i 

Wl  = 

64. 0  lb.' 

Fc  J  =  2.59  x  10% 

S3]  3 

+  1.138x10%  1 8, 


FDl-47.2  6 


-  5.7  X105  8  2 


for  |S3j  <  .12 

■ 

=(331. 2+5278. 0|83|) 

(83/|53|) 
for  |83[  i  .  12 


Figure  56.  Mass,  Spring  and  Damping  Force 
Equation  for  Final  Man  Model 


Mass  Spring  and  Damping  Constants  Used  In  Final  Non-  I 

Linear  Four- Degree-Of- Freedom  Model  | 


0 

1 

2 

3 

4 

Mode 

Legs/Feet 

Buttocks 

Visceral 

Spinal 

Head/Neck 

(Lower 

(Upper 

Pelvic 

Torso 

Mass) 

Mass) 

Weight  Distribution 

21 

64 

15 

48 

12 

(lb) 

(Total  Weight 

(Assumed 

=160  lb 

to  move  with 

139  lb 

feet) 

supported  by 

buttocks 

Spring  Constant 

Variable 

194 

Variable 

14,000 

k  (lb/ft) 

(eq.56) 

(eq.  31B) 

(eq.46) 

(eq.  42) 

Undamped 

natural  frequercy 

— 

47.5(1  g 

20.4 

65.0  for 

192.3 

,  4 

load) 

2-3  g  • 

variable 

load 

to  (rad/sec) 

(eq.54) 

(eq.31.A)  . 

(eq.45) 

(eq.4l) 

Spring  Force 

See  fig  55 

194  2 

See  fig  53 

14,000  4 

Fg  (lb) 

(eq.56) 

(eq.  31B) 

(eq.  46,47) 

(eq.42) 

Damping  Ratio  '  • 

.  25  (for 

0.5 

0.1  (for  3 

0.15 

c  (-) 

.  1=47.5) 

(eq.31.C) 

=65  rad/ 

(eq.39) 

(eq.57) 

sec) 

(eq.48) 

Damping 

Coefficient 

47.2 

9.5 

19.4 

10.75 

2K  (Ib/ft/sec) 

(eq.  57) 

(eq.31.  C) 

(eq.48) 

(eq.40) 

Damping  Force 

47.2 

9.5  2 

19.4  3 

10.75  . 

4 

Fd  (lb) 

(eq.57) 

(eq.  31.  C) 

(eq.48) 

(eq.  40) 

SECTION  V 


COMPARISON  BETWEEN  MODEL  OUTPUTS  AND  EXPERIMENT 


As  can  be  seen  from  figure  46,  linear  impedance  model  number  30  gives  good  agreement 
with  experiment  up  to  a  frequency  of  about  8  Hz.  Since  the  model  coefficients  were  ad¬ 
justed  with  this  end  in  view,  such  agreement  is  meaningless  by  itself,  of  course.  The 
lack  of  agreement  at  the  higher  frequencies  is  also  to  be  expected  because  of  the  import¬ 
ance  of  damping  at  the  higher  frequencies.  As  indicated  by  figure  10,  for  the  simple  case 
of  a  single-degree-of-freedom  cubic  model,  the  presence  of  nonlinear  damping  can  mar¬ 
kedly  change  the  apparent  impedance. 

If  impedance  model  number  30  (figure  46)  is  driven  by  a  sinusoidal  imput,  the  relative 
amplitudes  of  the  various  masses  are  as  indicated  in  figure  57.  These  results  are  very 
similar  to  the  experimental  measurements  of  Dieckmann  (ref.  13)  and  Latham  (ref.  14). 
Since  these  experiments  are  quite  different  to  impedance  measurements,  the  qualitative 
agreement  of  the  model  is  very  encouraging. 

Moving  now  to  the  final  model  configuration,  with  the  nonlinear'spine  and  buttock  springs, 
we  should  expect  similar  results  for  both  input  impedance  and  amplitude  ratio.  Unfor¬ 
tunately,  there  was  insufficient  time  available  to  carry  out  such  a  comparison  duringthe 
present  study. 

The  nonlinear  model  was  compared  with  an  AMRL  drop  test  tower  experiment,  however. 
As  indicated  in  figure  58,  the  input  acceleration  was  approximated  by  two  polynomials, 
and  the  resulting  force-time  history  is  given  in  figure  59.  Except  for  the  fact  that  the 
model  appears  to  be  slightly  under-damped,  the  agreement  with  experiment  is  considered 
to  be  excellent. 


ACCE  LEKATION  (in/see*) 


SECTION  VI 


SUGGESTED  FUTURE  WORK 


BUTTOCK  DYNAMICS 

The  stiffness  and  damping  of  the  buttock  mode  should  be  investigated  experimentally.  Some 
preliminary  work  in  this  area  has  been  carried  out  by  Phillips,  *  but  much  more  extensive 
testing  is  needed.  It  should  be  possible  to  establish  meaningful  methods  of  defining  buttock 
size  and  muscle  tone  and  to  correlate  the  dynamic  measurements  with  these  parameters. 
This  opens  the  way  to  determining  the  buttock  dynamic  characteristics  of  a  particular  air 
crew  member  by  simple  measurements  during  routine  physicals,  if  further  work  indicates 
that  buttock  dynamics  have  a  significant  influence  on  injury  potential. 

SPINAL  DYNAMICS  . 

It  should  be  relatively  easy  to  determine  spinal  stiffness  from  cadaver  material,  using 
complete  spines  from  subjects  of  different  ages,  and  there  seems  no  good  reason  why  this 
important  data  should  not  be  obtained  at  once.  Spinal  mode  damping  is  more  difficult,  since 
it  must  be  obtained  with  live  subjects.  At  the  moment,  and  excepting  the  work  of  Kazarian** 
with  monkeys,  even  the  mechanism  of  spinal  damping  is  unknown.  It  should  be  possible  to 
build  a  fixture  which  will  oscillate  the  upper  torso  with  respect  to  the  pelvic  girdle,  and  in 
this  way  obtain  damping  data  for  small  and  moderate  spine  deflections.  If  the  mechanism 
for  this  damping  can  then  be  discovered,  it  may  be  possible  to  extrapolate  the  damping 
data  throughout  the  range  of  operational  interest. 

GENERAL  RESEARCH  ON  DAMPING  MECHANISMS 

It  has  generally  been  assumed  that  damping  of  the  human  body  is  linear;  mostly,  one  sus¬ 
pects,  because  only  linear  damping  can  be  handled  analytically.  Kazarian's  work  with  mon¬ 
keys  implies  that  "orifice"  or  "hydraulic"  damping  may  be  significant  in  the  spinal  model. 
There  may  also  be  a, substantial  amount  of  quasi-coulomb  (constant  force)  damping  in  other 
parts  of  the  body,  however,  and  as  has  been  shown,  this  can  have  a  substantial  effect  on  the 
apparent  impedance  at  the  higher  frequencies.  . 

It  is  therefore  essential  to  obtain  a  better  understanding  of  the  mechanism  of  damping  in 
the  various  parts  of  the  human  body,  so  that  more  realistic  damping  relationships  can  be 
used  in  future  dynamic  models.  Some  new  possibilities  in  this  area  are  pointed  out  in  some 
earlier  work  done  by  Payne  and  Anthony.  *** 


*  See  footnote  on  p.  88. 

**  Verbal  communication  from  L.  Kazarian  of  AMRL. 

***  Unpublished  Wyle  Laboratory  working  paper  No.  140-1  (13  November  1969)  entitled 
"Notes  on  a  Simple  Model  of  a  Spinal  Vertebra"  by  P.  R.  Payne  and  A.  Anthony. 


RELATIVE  MOTION  OF  BODY  COMPONENTS 


Dieckmann  (ref.  13),  Latham  (ref.  14),  and  Woods  (ref.  15)  (among  others)  have  all  con¬ 
tributed  to  the  investigation  of  transmissibility,  and  the  relative  motion  of  body  components 
(head,  shoulder,  hips,  etc. )  under  steady  state  sinusoidal  excitation.  Unfortunately  their 
data  is  rather  fragmentary,  and  a  more  ambitious  program  needs  to  be  undertaken,  using 
at  least  ten  different  subjects,  with  very  carefully  designed  instrumentation.  It  might  be 
helpful  to  use  both  accelerometers  and  optical  techniques  at  the  same  time,  for  example. 

In  addition,  excitation  of  the  subject  should  be  applied  at  different  locations  and  bias,  as 
suggested  belsw. 

(a)  Excitation  at  zero  g  bias,  (with  the  subject  supported  horizontally  on  air  jets). 
Excitation  should  be  at  various  locations  (feet,  shoulders,  head,  for  the 
"standing"  position,  and  buttocks,  shoulders  and  head  for  the  "seated"  posi¬ 
tion). 

<b)  Similar  experiments  under  normal  1  g  bias. 

(c)  Similar  experiments  (excluding  head  excitation)  in  a  centrifuge,  to  obtain 
•  higher  bias  values. 

(d)  Experiments  where  "intermediate"  positions  (such  as  the  hips)  are  driven  . 
after  being  cast  in  plaster. 
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DISCUSSION  OF  AN  AMPLITUDE  RATIO  FUNCTION 
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Equating  to  zero  and  solving  for  p  at  df  =  q  we  obtain 

dt 

fy*?y  -*’* 

whence 


Substituting  back  into  the  equation  for  f  we  obtain,  after  some  manipulation 


Note  that  as 


1 


Note: 


O  (corresponding  to  zero  damping) 


/ 


Cx.~ 


(corresponds  to  jfZ,  ) 

(where  _f7res  is  resonant  frequency  for 

maximum  amplitude  ratio) 

« 

(corresponds  to  amplitude  ratio  -***.  c*a  ) 


Thus  for  low  damping  the  function  tends  to  behave  like  the  dynamic  magnifier  analyzed  in 
Appendix  II,  near  the  point  of  resonance. 

A  typical  plot  of  the  function  f  is  given  in  figure  1. 1.  The  maximum  value  of  f,  the  corres¬ 
ponding  value  of  the  frer  iency  ratio  parameter  p,  is  plotted  in  figure  1.2. 
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APPENDIX  III 


NOTE  ON  COMPUTER  TECHNIQUES  EMPLOYED  IN 
LUMPED  PARAMETER  SYSTEM  ANALYSIS 

,  /  ■ 

•  *  . 

During  the  analysis  described  in  this  rep^«  a  number  of  computational  techniques  were 
employed. 

LINEAR  MODELS 


The  simplest  method  of  treating  the  linear  models  was  found  to  be  by  use  of  a  simple  elec¬ 
tric  analog  as  shown  in  figure  III.  1.  The  one-to-one  correspondence  between  electrical 
and  mechanical  impedance  and  between,  respectively,  inductance  and  mass,  capacitance 
and  elasticity,  and  resistance  and  damping  coefficient  coupled  with  the  existence  of  electric 
circuit  analysis  routines  made  the  analysis  of  these  cases  extremely  easy.  This  method 
was  used  for  developing  the  frequency  impedance  curves  shown  in  figures  22  to  51.  This 
method  is  currently  not  flexible  enough"  to  accept  arbitrary  transient  inputs. 


NONLINEAR  SPRING  MODELS 


The  cases  illustrated  in  figure!  53  through  59  in  which  the  buttock  and  spinal  spring  forces 
were  represented  by  nonlinear  discontinuous  functions  could  not  be  treated  by  the  simple 
electric  analog  method  used  for  the  linear  model. 


In  the  nonlinear  cases  a  digital  simulation  of  a  75  block  analog  program  was  used.  This 
program  was  sufficiently  flexible  to  allow  the  introduction  of  a  transient  input  acceleration 
such  as  that  shown  in  figure  58. 


MODELS  WITH  NONLINEAR  DAMPING  AND  SPRINGING 

The  case  (illustrated  in  figure  10)  in  which  nonlinear  damping  and  springing  was  used  was 
programmed  in  a  conventional  digital  program.  A  Runge-Kutta  step-by-step  integration 
routine  was  used  and  sufficient  cycles  were  calculated  to  ensure  a  steady-state  response 
had  been  reached.  This  digital  program  could  readily  be  adapted  to  the  constant  friction 
force  case  or  to  transient  inputs  although  the  latter  were  not  actually  used. 
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